SupeRNAlign: a new tool for flexible superposition of homologous RNA structures and inference of accurate structure-based sequence alignments

RNA has been found to play an ever-increasing role in a variety of biological processes. The function of most non-coding RNA molecules depends on their structure. Comparing and classifying macromolecular 3D structures is of crucial importance for structure-based function inference and it is used in the characterization of functional motifs and in structure prediction by comparative modeling. However, compared to the numerous methods for protein structure superposition, there are few tools dedicated to the superimposing of RNA 3D structures. Here, we present SupeRNAlign (v1.3.1), a new method for flexible superposition of RNA 3D structures, and SupeRNAlign-Coffee—a workflow that combines SupeRNAlign with T-Coffee for inferring structure-based sequence alignments. The methods have been benchmarked with eight other methods for RNA structural superposition and alignment. The benchmark included 151 structures from 32 RNA families (with a total of 1734 pairwise superpositions). The accuracy of superpositions was assessed by comparing structure-based sequence alignments to the reference alignments from the Rfam database. SupeRNAlign and SupeRNAlign-Coffee achieved significantly higher scores than most of the benchmarked methods: SupeRNAlign generated the most accurate sequence alignments among the structure superposition methods, and SupeRNAlign-Coffee performed best among the sequence alignment methods

Genome-wide survey of codons under diversifying selection in a highly recombining bacterial species, Helicobacter pylori

Selection has been a central issue in biology in eukaryotes as well as prokaryotes. Inference of selection in recombining bacterial species, compared with clonal ones, has been a challenge. It is not known how codons under diversifying selection are distributed along the chromosome or among functional categories or how frequently such codons are subject to mutual homologous recombination. Here, we explored these questions by analysing genes present in >90% among 29 genomes of Helicobacter pylori, one of the bacterial species with the highest mutation and recombination rates. By a method for recombining sequences, we identified codons under diversifying selection (dN/dS > 1), which were widely distributed and accounted for ∼0.2% of all the codons of the genome. The codons were enriched in genes of host interaction/cell surface and genome maintenance (DNA replication, recombination, repair, and restriction modification system). The encoded amino acid residues were sometimes found adjacent to critical catalytic/binding residues in protein structures. Furthermore, by estimating the intensity of homologous recombination at a single nucleotide level, we found that these codons appear to be more frequently subject to recombination. We expect that the present study provides a new approach to population genomics of selection in recombining prokaryotes

Automated functional classification of experimental and predicted protein structures

BACKGROUND: Proteins that are similar in sequence or structure may perform different functions in nature. In such cases, function cannot be inferred from sequence or structural similarity. RESULTS: We analyzed experimental structures belonging to the Structural Classification of Proteins (SCOP) database and showed that about half of them belong to multi-functional fold families for which protein similarity alone is not adequate to assign function. We also analyzed predicted structures from the LiveBench and the PDB-CAFASP experiments and showed that accurate homology-based functional assignments cannot be achieved approximately one third of the time, when the protein is a member of a multi-functional fold family. We then conducted extended performance evaluation and comparisons on both experimental and predicted structures using our Functional Signatures from Structural Alignments (FSSA) algorithm that we previously developed to handle the problem of classifying proteins belonging to multi-functional fold families. CONCLUSION: The results indicate that the FSSA algorithm has better accuracy when compared to homology-based approaches for functional classification of both experimental and predicted protein structures, in part due to its use of local, as opposed to global, information for classifying function. The FSSA algorithm has also been implemented as a webserver and is available at

Evaluation of 3D-Jury on CASP7 models

<p>Abstract</p> <p>Background</p> <p>3D-Jury, the structure prediction consensus method publicly available in the Meta Server <url>http://meta.bioinfo.pl/</url>, was evaluated using models gathered in the 7^{<it>th </it>}round of the Critical Assessment of Techniques for Protein Structure Prediction (CASP7). 3D-Jury is an automated expert process that generates protein structure meta-predictions from sets of models obtained from partner servers.</p> <p>Results</p> <p>The performance of 3D-Jury was analysed for three aspects. First, we examined the correlation between the 3D-Jury score and a model quality measure: the number of correctly predicted residues. The 3D-Jury score was shown to correlate significantly with the number of correctly predicted residues, the correlation is good enough to be used for prediction. 3D-Jury was also found to improve upon the competing servers' choice of the best structure model in most cases. The value of the 3D-Jury score as a generic reliability measure was also examined. We found that the 3D-Jury score separates bad models from good models better than the reliability score of the original server in 27 cases and falls short of it in only 5 cases out of a total of 38. We report the release of a new Meta Server feature: instant 3D-Jury scoring of uploaded user models.</p> <p>Conclusion</p> <p>The 3D-Jury score continues to be a good indicator of structural model quality. It also provides a generic reliability score, especially important for models that were not assigned such by the original server. Individual structure modellers can also benefit from the 3D-Jury scoring system by testing their models in the new instant scoring feature <url>http://meta.bioinfo.pl/compare_your_model_example.pl</url> available in the Meta Server.</p

Circular Permutation in Proteins

This is a ‘‘Topic Page’ ’ article for PLoS Computational Biology. Circular permutation describes a type of relationship between proteins, whereby the proteins have a changed order of amino acids in their protein sequence, such that the sequence of the first portion of one protein (adjacent to the N-terminus) is related to that of the second portion of the other protein (near its C-terminus), and vice versa (see Figure 1). This is directly analogous to the mathematical notion of a cyclic permutation over the set of residues in a protein. Circular permutation can be the result of evolutionary events, post-translational modifications, or artificially engineered mutations. The result is a protein structure with different connectivity, but overall similar three-dimensional (3D) shape. The homology between portions of the proteins can be established by observing similar sequences between N- and C-terminal portions of the tw

A Freeze Frame View of Vesicular Stomatitis Virus Transcription Defines a Minimal Length of RNA for 5′ Processing

The RNA synthesis machinery of vesicular stomatitis virus (VSV) comprises the genomic RNA encapsidated by the viral nucleocapsid protein (N) and associated with the RNA dependent RNA polymerase, the viral components of which are a large protein (L) and an accessory phosphoprotein (P). The 241 kDa L protein contains all the enzymatic activities necessary for synthesis of the viral mRNAs, including capping, cap methylation and polyadenylation. Those RNA processing reactions are intimately coordinated with nucleotide polymerization such that failure to cap results in termination of transcription and failure to methylate can result in hyper polyadenylation. The mRNA processing reactions thus serve as a critical check point in viral RNA synthesis which may control the synthesis of incorrectly modified RNAs. Here, we report the length at which viral transcripts first gain access to the capping machinery during synthesis. By reconstitution of transcription in vitro with highly purified recombinant polymerase and engineered templates in which we omitted sites for incorporation of UTP, we found that transcripts that were 30-nucleotides in length were uncapped, whereas those that were 31-nucleotides in length contained a cap structure. The minimal RNA length required for mRNA cap addition was also sufficient for methylation since the 31-nucleotide long transcripts were methylated at both ribose-2′-O and guanine-N-7 positions. This work provides insights into the spatial relationship between the active sites for the RNA dependent RNA polymerase and polyribonucleotidyltransferase responsible for capping of the viral RNA. We combine the present findings with our recently described electron microscopic structure of the VSV polymerase and propose a model of how the spatial arrangement of the capping activities of L may influence nucleotide polymerization

Distant homologs of anti-apoptotic factor HAX1 encode parvalbumin-like calcium binding proteins

<p>Abstract</p> <p>Background</p> <p>Apoptosis is a highly ordered and orchestrated multiphase process controlled by the numerous cellular and extra-cellular signals, which executes the programmed cell death <it>via </it>release of cytochrome c alterations in calcium signaling, caspase-dependent limited proteolysis and DNA fragmentation. Besides the general modifiers of apoptosis, several tissue-specific regulators of this process were identified including HAX1 (HS-1 associated protein X-1) - an anti-apoptotic factor active in myeloid cells. Although HAX1 was the subject of various experimental studies, the mechanisms of its action and a functional link connected with the regulation of apoptosis still remains highly speculative.</p> <p>Findings</p> <p>Here we provide the data which suggests that HAX1 may act as a regulator or as a sensor of calcium. On the basis of iterative similarity searches, we identified a set of distant homologs of HAX1 in insects. The applied fold recognition protocol gives us strong evidence that the distant insects' homologs of HAX1 are novel parvalbumin-like calcium binding proteins. Although the whole three EF-hands fold is not preserved in vertebrate our analysis suggests that there is an existence of a potential single EF-hand calcium binding site in HAX1. The molecular mechanism of its action remains to be identified, but the risen hypothesis easily translates into previously reported lines of various data on the HAX1 biology as well as, provides us a direct link to the regulation of apoptosis. Moreover, we also report that other family of myeloid specific apoptosis regulators - myeloid leukemia factors (MLF1, MLF2) share the homologous C-terminal domain and taxonomic distribution with HAX1.</p> <p>Conclusions</p> <p>Performed structural and active sites analyses gave new insights into mechanisms of HAX1 and MLF families in apoptosis process and suggested possible role of HAX1 in calcium-binding, still the analyses require further experimental verification.</p

Editing of hnRNP K protein mRNA in colorectal adenocarcinoma and surrounding mucosa

The heterogeneous nuclear ribonucleoprotein K (hnRNP K) protein is an RNA-binding protein involved in many processes that compose gene expression. K protein is upregulated in the malignant processes and has been shown to modulate the expression of genes involved in mitogenic responses and tumorigenesis. To explore the possibility that there are alternative isoforms of K protein expressed in colon cancer, we amplified and sequenced K protein mRNA that was isolated from colorectal cancers as well as from normal tissues surrounding the tumours. Sequencing revealed a single G-to-A base substitution at position 274 that was found in tumours and surrounding mucosa, but not in individuals that had no colorectal tumour. This substitution most likely reflects an RNA editing event because it was not found in the corresponding genomic DNAs. Sequencing of RNA from normal colonic mucosa of patients with prior resection of colorectal cancer revealed only the wild-type K protein transcript, indicating that G274A isoform is tumour related. To our knowledge, this is the first example of an RNA editing event in cancer and its surrounding tissue, a finding that may offer a new diagnostic and treatment marker